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We exhaust all exact self-consistent solutions of complex-valued fermionic condensates in the 1 + 1 dimen- 
sional Bogoliubov-de Gennes and chiral Gross-Neveu systems under uniform boundary conditions. We obtain 
n complex (twisted) kinks, or grey solitons, with 2n parameters corresponding to their positions and phase 
shifts. Each soliton can be placed at an arbitrary position while the self-consistency requires its phase shift to 
be quantized by n/N for N flavors. 
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INTRODUCTION 

The search for inhomogeneous self-consistent fermionic 
condensates including states such as the Fulde-Ferrell (FF) 
[1] and Larkin-Ovchinnikov (LO) [2] states having phase and 
amplitude modulations, respectively, in superconductors has 
attracted considerable attentions for more than half a century 
after theoretical predictions were made about their existence. 
While amplitude modulations are well studied in conduct- 
ing polymers [3-7], the FFLO state is theoretically shown to 
be a ground state of superconductors under a magnetic field 
[8]. Recently, the FFLO state has also been discussed in the 
context of ultracold atomic Fermi gas [9, 10]. In general, 
inhomogeneous self-consistent fermionic condensates with a 
gap function and fermionic excitations can be treated simul- 
taneously using the Bogoliubov-de Gennes (BdG) and gap 
equations [11]. The gap functions are real and complex for 
conducting polymers [12] and superconductors, respectively. 
In the quantum field theory, these systems correspond to the 
Gross-Neveu (GN) model [13] and the Nambu-Jona-Lasinio 
(or chiral GN) model [14], which were proposed as models 
of dynamical chiral symmetry breaking in 1+1 or 2+1 dimen- 
sions. Therefore, BdG and (chiral) GN systems have been 
studied and developed together from the viewpoint of both 
condensed matter physics and high energy physics (see [15] 
as a review). For instance, fermion number fractionization is 
one of the topics that has been studied from this viewpoint 
[16]. Recently, it has been shown that the solutions in 1 + 1 
dimensions can be promoted to 3+1 dimensions [17], thereby 
leading to extensive study of the modulated phases of these 
systems in terms of quantum chromodynamics (QCD) [18]. 

Inhomogeneous self-consistent solutions are often studied 
numerically because analytic solutions are generally difficult 
to obtain. However, several analytic solutions are available 
in the case of the real-valued condensates in 1 + 1 dimen- 
sions, which describe the conducting polymers and the real 
GN model. Under uniform boundary conditions at spatial in- 
finities, a real kink was constructed by Dashen et. al. [19] by 
using the inverse scattering method, and later, it was recon- 



structed in polyacetylene [20] in the continuum limit of the 
lattice model [21]. Subsequently, a bound state of a kink and 
an anti-kink, which is called a polaron, was constructed in 
polyacetylene [22, 23], for which achieving self-consistency 
in the system requires the distance between the kink and anti- 
kink to be fixed. Furthermore, three kinks (kink and polaron 
placed at arbitrary positions) [24, 25] and more general so- 
lutions [26] were obtained. The attractive interaction between 
two polarons was also investigated [27]. For a periodic bound- 
ary condition, the existence of real kink crystals (the LO state) 
has been known for a long time [3-7]. 

On the other hand, when compared with real conden- 
sates, only a few self-consistent solutions have thus far been 
obtained for complex condensates, such as a complex (or 
twisted) kink or a grey soliton, [28] and their crystals [29, 30]. 
In these complex-valued crystals, both the amplitude and 
phase are modulated (the FFLO state), and this modulated 
phase has important applications in both superconductors 
and QCD, such as in the phase diagram of the chiral GN 
model [31]. An attempt to construct more general solutions 
was made [32, 33] by using a technique of integrable sys- 
tems known as the nonlinear Schrodinger or Ablowitz-Kaup- 
Newell-Segur hierarchy [34]. 

In this Letter, we exhaust all exact self-consistent solutions 
of complex condensates under uniform boundary conditions, 
and we find that they describe multiple twisted kinks. Unlike 
polarons in real condensates, where the distance between the 
kink and anti-kink is fixed, the situation is drastically simpli- 
fied in our multiple twisted-kink solutions; we determine the 
filling rate of fermions for bound states of each kink, and we 
find that each kink can be placed at any position and has any 
phase shift quantized by n/N with the number of flavors N. 



RESULT FROM THE INVERSE SCATTERING THEORY 

First, we briefly summarize the mathematical expressions 
of the n-soliton solution and its eigenstates of the self- 
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FIG. 1. Scattering and bound states in s-plane. Scattering states with 
positive (negative) energy exist on the real and positive (negative) 
axis. Bound states exist on the unit circle, and s and s* represent the 
same bound state. 
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FIG. 2. Example of 3-kink solution. Here the parameters are S\ 
eA m , s 2 = ei" 1 , s 3 = ei* 1 , X\ = -10, x 2 = 10, and x 3 = 0. 



defocusing Zakharov-Shabat eigenvalue problem [35] 



-id x A(x) 
AW* id x 



(i) e 



obtained by the inverse scattering method [36]. 

Uniformizing variable. Let us assume that the gap function 
obeys the following asymptotically uniform boundary condi- 
tion: 



|A(x)| -> m (> 0), 



(2) 



Subsequently, we parametrize the energy and wavenumber of 
quasiparticles using the uniformizing variable s defined by 



e(j) = -(s + s- 1 ), 



)■ 



(3) 



k 2 + m 2 



We can easily verify that the dispersion relation e 2 
holds for an arbitrary complex number s. Eigenstates corre- 
sponding to s on the real axis are scattering states, while those 
on the unit circle are bound states. Since s and s* on the unit 
circle correspond to the same bound state, it is sufficient to 
consider the unit circle in the upper half-plane when we count 
the number of bound states (Figure 1). The ranges of the vari- 
able are as follows 



-oo < s < <-> e<0&-oo<fc<oo, 
< < oo <-> e > & — oo<fc<oo, 



(4) 



the 



and therefore when we rewrite the integral of k using s, 
interval becomes (-oo, 0) for e < and (0, +oo) for e > 0. 

n-soliton solution. Let us consider the gap function A(x) 
which has n bound states and acts as a reflectionless potential 
for scattering states, i.e., the n-soliton solution. By writing the 



s values of the bound states as s; 



(J = !,...,«) with 



< Qj < 7T, the eigenenergy and the complex wavenumber 
can be rewritten as 



Kj :- -ik(sj) = m sin#/ 



£j := c(sj) = mcosOj, (5) 



respectively. According to the inverse scattering theory, 
i , . . . , 9„ are all different from each other and there is no de- 
generacy. We further introduce the following variables: 



ej(x) = e' 



(7 = l ( ...,n). 



(6) 



Here, the real constant xj roughly represents the position of 
the j-th soliton. (It is to be noted that sj values also describe 
the positions of the solitons.) Furthermore, we define the func- 
tions f\ (x), . . . , f n (x) as solutions of the following linear equa- 
tion: 
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= 0. (7) 



Here, the argument x is abbreviated. 

By using the above notations, the n-soliton solution can be 
expressed as 



A(x) = m + 2i 2] Sj l ej(x)fj(x). 



(8) 



./=! 



The complex-valued n-soliton solution has 2n parameters 
si , . . . , s n , x\ , . . . , x„, and this number of parameters is exactly 
twice that of the real-valued soliton solution. This A(x) has 
the following asymptotic form: 



A« 



(x 



3), 



-2i(e l +B 2 +-+e„) 



(x — > +00). 



(9) 



If the solitons are sufficiently separated from each other, the 
phase shift brought about by the 7-th soliton is sf = e~ 2l0 K 
Figure 2 shows an example of the 3-soliton solution. 

The reduction to the real-valued soliton solution is obtained 
as follows. When the number of solitons is even (n = 2n'), the 
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relations 



S2J-1 



X2j-i =x 2 j (j = !,...,«') (10) 



yield real-valued solutions. When the number of solitons is 
odd (« = 2n' + 1), we need to consider the term si n '+i = & m ^ 2 
in addition to Eq. (10) while X2 n >+i remains arbitrary. By 
this reduction, we obtain f 2 j-i(x) = f2j(x)* and f2 n '+i(x) = 
f2n'+\{x)*, and the imaginary part of Eq. (8) vanishes. 
The bound state with s = sj (<-> e = m cos dj) is given by 



Uj(x)\ ( fj(x) 
Vj( X )J \sjfj(xy)> 



(j= !,...,«). 



(11) 



We can show that this state is already normalized, i.e., 
J dx(\uj\ 2 + \ Vj \ 2 ) = 1 holds. 

Finally, let s be real. Then, the scattering states are given 

by 



lu(x, s)\ _ ^ 
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zi ^ 
m £—1 



7=1 J 



e A x ) I fj(x) 
S \Sjfj(x)' 



(12) 



which are obviously reflectionless as observed from the ex- 
pression. The amplitudes of these solutions at x = ±oo are 



[w(+oo, s)\ + [v(±oo, s)\ - 1 + s 



(13) 



A detailed derivation of these expressions for eigenstates will 
be published elsewhere. 



GAP EQUATION 

The one-dimensional BdG system with the Andreev ap- 
proximation consists of the BdG equation for right movers 
(BdG R ) 



id x A(x)\ ( U R \ I Mr 



\A(xy id x J\v R ) \v R 
the BdG equation for left movers (BdGO 



id. v A(x)\ ( u h \ _ tu h 
A(x)* -idJW ~ e \v L 



and the gap equation as a self-consistent condition 

A(x) 



2^ (urv* r + u L v* L ) . 

occupied states 



(14) 



(15) 



(16) 



The following relation exists between the solutions of the 
right and left movers: 



(e, u(x), v(x)) is a solution of BdGR. 
<-> (-€*, -v(x)*, u(x)*) is a solution of BdGL- 



(17) 



Thus, we can rewrite all quasiparticle wavefunctions of the 
left movers using those of the right movers. 

From this point onwards, we consider the BdG system 
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FIG. 3. Diagram of the occupation states considered in this Letter. In 
this example figure, the number of flavors is N = 6 and the number 
of solitons is n = 3. The filling rates defined by Eq. (19) are given 
by Vi = -1/6, v 2 = 1/3, and v 3 = 1/2. 



with N internal degrees of freedom, or equivalently, the chi- 
ral GN model with N flavors. In the light of examining low- 
energy excited states of condensed matter systems, we con- 
sider the configurations in which all the negative-energy scat- 
tering states are filled by fermions and positive-energy states 
are completely vacant. As for bound states, we label the bound 
states of BdGR as (m ^r, v^r) (J = l,...,ri), and we also label 
the corresponding bound states of BdGL with the energy of 
the opposite sign related by (17) as («j,l>v,;l)- These states 
are assumed to be filled partially, and we write the occupation 
number as Nj R and Njl, as schematically shown in Figure 3. 
The gap equation subsequently becomes 

A(jc) 
g 

= Nu r v* r + 2 Nu L vl + 2 N M u jtR v* jtR + Y A',, r). 



£<(> 



£<(> 



N( «rv r - «R v r) + X ( ^' R " N J0 U J» V U- (18) 
l-o l>o I bs - 



Here, the notation s.s. (b.s.) denotes the scattering (bound) 
states, and the relation (17) is used to obtain the second line 
in the above equation. Defining the filling rate by 



N jR -N jL 



■l<vj<l, (j =!,...,«), (19) 



1 ' N 

the above equation can be rewritten as follows: 

A(jc) 



=- = 2 MrV R ~ MrV R + Tj y j u J- RV U (20) 



S.S. S.S. K s 

e<0 e>0 



with g := Ng. It is to be noted that the sum of positive-energy 
scattering states in Eq. (20) has a negative sign because of the 
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relation (17), and it is equivalent to the stationary condition of 
the action in the GN model, as given in Ref. [28]. Thus, we 
can again confirm the equivalence of the problems between 
the BdG and GN systems. 

Henceforth, we always use the quantities of the BdGR 
system, and we omit the subscript R. Considering the limit 
L — > oo, where L denotes the system size, we replace the 
sum of the scattering states of the gap equation (20) by the 
corresponding integral. After renormalization of the coupling 
constant 



lip dk 

g 2n Xoo V£ 2 + m 



(21) 



and subtracting the logarithmically divergent terms from both 
sides, we obtain the following expression: 



= Yj v i u A x ) v A x Y 



i r°° dk I A(x) 



) u k {x)v k {xT 
-(sgne)- ij— 2| , 



(22) 



where we have written the scattering states with the wavenum- 
ber k as {uk{x), Vk(x)), and their amplitudes at infinity as 
("yt.co, Vjt oc). (We note that (u^(x), v k (x)) for a positive en- 
ergy and that for a negative energy are different from each 
other, though we use the same notation.) It is convenient to 
rewrite the above integral in terms of the uniformizing vari- 
able s introduced in Eq. (3). Using the relation (4) and 
\s\(l + s~ 2 ), we obtain 



V/t 2 + m 2 



= X v i u M) v A x Y 

b.s. 

r° pi ds {mil + s- 2 )u(x, s)v(x, sf 

J-oo Jo 



A(x) 

2tt I 2(|k(oo, s)\ 2 + | v(oo, s )\ 2 ) 2s 



(23) 



Here, we have written the scattering states labeled by s as 

(u(x, s), v(x, s)). 



SELF-CONSISTENT CONDITION FOR COMPLEX 
n-SOLITON SOLUTION 

We first present our main result in the following theorem, 
and then provide the proof. 

Theorem. Let A(x) be an w-soliton solution given by Eq. (8). 
The gap equation (23) holds if and only if the filling rate Vj 
satisfies 



26 j - n 



(24) 



Here, we remark on certain aspects of this theorem: 
1 . This theorem provides all self-consistent solutions under 
the uniform boundary condition (2), because A(x) needs to be 
a reflectionless potential in order for the gap equation to hold 
[19, 26, 28], and n-soliton solutions cover all reflectionless 
potentials with n bound states. 



2. The filling rate v, for the y'-th bound state only depends 
on the phase shift of the y'-th soliton, and it is not affected by 
other soliton parameters. Thus, the self-consistent condition 
is decoupled for each bound state (or each soliton). 

3. The position of the soliton Xj {j — 1, . . . , n) is arbitrary and 
is not related to the self-consistency. This contrasts with the 
case of real-valued condensates. Because they must be real, 
the distance between two solitons must be fixed to a specific 
value, such as in the case of the polarons in polyacetylene [22- 
24] and the topologically trivial soliton in the GN model [25]. 

4. For the A^-flavor system, the possible values of the filling 
rate are given by v = , ^-p . Correspondingly, 
the possible phase shift of each soliton is also discretized. For 
example, only the trivial value v = is allowed for N = 1, 
which corresponds to the real kink 29 = ji. When N — 2, 
the values v = -A, 0, i are allowed, which correspond to 
6 = |, |, *f. The cases N — 1 and 2 correspond to .s-wave 
superconductors and polyacetylene, respectively. On the other 
hand, any soliton solution can be self-consistent when N = oo. 

Proof. Upon substituting Eqs. (8), (12), and (13) into the 
integrand of the gap equation (23), the terms 1 + s~ 2 and 
\u(±oo, s)| 2 + |v(±oo, s)\ 2 in the first term in the bracket undergo 
cancellation, thereby yielding 

A(x) 



—u(x, s)v(x, s)* 
2 2s 



sinfy 4i y~< e jfj e ifi sin#/ 



-?V f j V jJJ '■ 



(25) 



j jJ 
From this equation and Eq. (7) expressed as 



si \s - Sj[ 



= 0, 



C/=l,...,n), 



we obtain the integrand of the gap equation (23) as 
A(x) n ^-i _ lx2 sinfy 



—u(x, s)v(x, s) - — — =-2^ s j fj 



For the elementary formula 
ds 1 



Ji 



, ,12 



tan 



sinf 



sin 0j 



(26) 



(27) 



(28) 



we can perform the integration to obtain 



sY- 



MxY 

2s i 



26i - n 



(29) 



Recalling that the bound states are given by Eq. (1 1), we fi- 
nally obtain 



(R.H.S. of Eq. (23)) = £ sj 1 f 2 (v, - ^-^J . 



(30) 



Since the functions fi(x) 2 , . . ., f„(x) 2 are linearly independent 
of each other, the theorem holds. □ 
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SUMMARY 

In summary, we have constructed all the exact self- 
consistent solutions of complex condensates under uniform 
boundary conditions. Our multiple n-twisted kink solution 
contains 2n parameters, and each kink has one bound state. 
Each kink can be placed at any position, while the self- 
consistency of the system requires the phase shift of each kink 
to be quantized by n/N with the number of flavors N. Our 
solution describes multiple grey solitons in ultracold atomic 
fermion gases, and our predictions require experimental veri- 
fication. The dynamics and scattering of these solitons should 
be studied as a future research topic. Further research also 
needs to be conducted on the construction of self-consistent 
solutions under non-uniform boundary conditions. 
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